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Abstract 

We demonstrate that hydrodynamic reductions of dispersionless integrable systems in 
2 + 1 dimensions, such as the dispersionless Kadomtsev-Petviashvili (dKP) and dispersionless 
Toda lattice (dTl) equations, can be deformed into reductions of the corresponding dispersive 
counterparts. Modulo the Miura group, such deformations are unique. The requirement that 
any hydrodynamic reduction possesses a deformation of this kind imposes strong constraints 

SlJ I on the structure of dispersive terms, suggesting an alternative approach to the integrability 

^y. • in 2 + 1 dimensions. 
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1 Introduction 

The Kadomtsev-Petviashvili (KP) equation, 

arises in mathematical physics as a two-dimensional generalization of the KdV equation. Al- 
though its integrability aspects have been thoroughly investigated in the literature, we believe 
that one important property has been overlooked, namely, that Eq. ([T|) can be decoupled into 
a pair of consistent (1 -|- 1) -dimensional equations in a continuum of ways. These decouplings 
can be obtained as deformations of hydrodynamic reductions of its dispersionless limit, known 
as the dKP equation, 

{Ut - UUx)x = Uyy (2) 

We point out that Eq. ^, also known as the Khokhlov-Zabolotskaya equation [3_1J, is of interest 
in its own, arising in non-linear acoustics, gas dynamics and differential geometry. A key property 
of the dKP equation is the existence of n-phase solutions of the form 

u = u{R^,...,R''), (3) 

where the 'phases' R''{x,y,t) are governed by a pair of commuting hydrodynamic type systems 

Rl = fi\R)Ri, Rl = X'{R)Ri. (4) 

Here A* = (/u*)^ + u, while //* and u satisfy the so-called Gibbons-Tsarev equations, 

a,.. = ^, 3.3,u = 2j^^, (5) 

^ 7^ J) ^j = d/dK^, which were first derived in |151ll6j in the context of hydrodynamic reductions 
of Benney's moment equations. Thus, Eqs. (|1|) can be viewed as a decomposition of the (2 + 1)- 
dimensional Eq. ([2]) into a pair of (1 + l)-dimensional hydrodynamic type systems. We will refer 
to Eqs. ([!]) as hydrodynamic reductions of dKP. These reductions have been extensively studied 
in the literature, see e.g. |141 [HI [22] and references therein. In particular, in the one-component 
case Eqs. ([5]) become vacuous, and without any loss of generality one can set u{R) = R where 
R{x,y,t) satisfies a pair of Hopf-type equations 

Ry = fiRx, Rt = ifJ' + R)Rx', (6) 

here n{R) is an arbitrary function. We recall that the general solution of Eqs. ([6|) is given by 
the implicit formula f{R) = x + ^y + {p? + R)t, which implies that the level surfaces i?=const 
are planes. Solutions of this type are known as planar simple waves. 

Our main observation is that all hydrodynamic reductions (jl]) can be deformed into reduc- 
tions of the full KP equation by adding appropriate dispersive terms which are polynomial in the 
x-derivatives of -R*. Up to Miura-type transformations, such deformations are unique. Moreover, 
the calculation of dispersive corrections is an entirely algebraic procedure which does not require 
solving differential equations. In the one-component case one obtains the following deformation 
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ofEqs. ®: 

y — f^^x 

+ ^ ((2^^' + 1)R,, + (^/i" - fiif^'f + i^rmRl), + 0{e% 

see Sect. 2.1 for more details. Notice that the relation u = R, which now solves the full KP 
equation, remains undeformed: this can always be assumed modulo the Miura group. Eqs. ([7]) 
can be viewed as a decomposition of KP into a pair of commuting (1 + l)-dimensional equations 
parametrized by an arbitrary function of one variable. In general, these equations constitute in- 
finite series in e^ which can terminate only in exceptional cases (in the one-component situation 
this happens only when /U=const, in which case KP reduces to KdV). One can show that defor- 
mations ([7]) are nontrivial, that is, not reducible to Eqs. ([6]) by a Miura-type transformation. 
Formal expansions of the type ([7]) have appeared in [1] , and were thoroughly investigated in a 
series of publications [SJ El [TJ [2T| [28] in the context of 2D topological field theory. We would 
like to formulate the following conjecture: 

For any integrable soliton system in 2 + 1 dimensions, all hydrodynaniic reductions of its diper- 
sionless limit can be deformed into reductions of the original system. 

Thus, any (2 + l)-dimensional integrable soliton equation can be decoupled into a pair of com- 
patible (1 + l)-dimensional equations in an infinity of ways. Deformations of two-component 
reductions of dKP are discussed in Sect. 2.2 and 2.3. Another example supporting our conjecture 
is provided by the Toda lattice system, 

euy = u{w{x)-w{x-e)), 
ewt = u{x + e) — u{x), 

whose dispersionless limit (dTl equations) assumes the form 

Uy = UWx, Wt = Ux- 

In Sect. 3 we provide explicit formulae for dispersive deformations of one-component reductions 
of dTl. 

The requirement that all hydrodynamic reductions of a (2 + l)-dimensional dispersionless 
integrable system can be deformed into reductions of the corresponding dispersive equation 
imposes strong constraints on the structure of dispersive terms. As an illustration, let us consider 
the generalized KP equation of the form 

{ut - uux + e{AiUxx + Mul) + e^{BiUxxx + B2UxUxx + B3ul))^ = Uyy, (9) 

where Ai and Bi are certain functions of u. Notice that this equation has the same dispersionless 
limit as the KP equation ([T]). As demonstrated in Sect. 4.1, the requirement that all one- 
component reductions of dKP can be deformed into reductions of ([9|) readily implies Ai = A2 = 
B2 = B-^ = and i?i=const, moreover, to establish this one only needs to perform calculations up 
to the order e^. Thus, our procedure uniquely reconstructs the KP equation. Further examples 



of this type include BKP/CKP and the 'universal' equation (Sect. 4.2 and 4.3, respectively). 
This naturally leads to the program of classification of (2 + l)-dimensional integrable soliton 
equations which can be summarized as follows: 

(a) Classify (2 + l)-dimensional dispersionless integrable systems within various particularly 
interesting classes. A number of results in this direction are already available, see e.g. [21 [HI 
[HI [ini [m [121 [13], etc. We recall that the integrability of a (2 + l)-dimensional dispersionless 
system is understood as the existence, for any n, of an infinity of n-component hydrodynamic 
reductions parametrized by n arbitrary functions of one variable [S] . 

(b) Reconstruct possible dispersive terms from the requirement that all hydrodynamic reductions 
of the dispersionless system can be deformed into reductions of the corresponding dispersive 
equation. We conjecture that any (2 + l)-dimensional dispersionless integrable system can be 
deformed in this way (possibly, non-uniquely) . 

This scheme can be viewed as an alternative to the classical approach to the integrability 

2 

in 2 + 1 dimensions, which starts with a linear dispersive part, say, utx — f^'^xxxx = Uyy, and 
reconstructs the allowed nonlinearity 



2 Deformations of the dKP reductions 

We will consider separately deformations of one- and two-component reductions of dKP. 

2.1 One-component reductions 

Let us first rewrite both KP and dKP equations as two-component systems, 

£2 

Ut - UUx - —UxXX = Wy, Uy = Wx (lO) 

and 

Ut - UUx = Wy, Uy = Wx, (11) 

respectively; such representation simplifies the calculations. One-component (one-phase) reduc- 
tions of dKP are given by the formula 

u = R, w = w{R), (12) 

where w' = fi and R{x, y, t) solves Eqs. ([6|). Dispersive deformations of Eqs. ([T2]) . ([6|) are sought 
in the form 



u =R, 

w =w{R) 

+ €\biRxx + b2Rl) 

+ e {diRix + d2RxxxRx + d^R^^ + d^RxxRx + d^Rx) + 0{e ), 



(13) 



where the coefficients bi and di are certain functions of R, and d^R = Rnx- Notice that the 
relation u = R remains undeformed: this can always be achieved modulo the Miura group. 
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Similarly, the deformed version of Eqs. © is 
Ry =iiRx 

+ €^{ipiRxxx + ^2RxxRx + f-iRl) 

+ e {PlRbx + P2RaxRx + PsRxxxRxx + P^RxxRx + P^RxxxRx + PqRxxRx + PlRx) 

+ 0{e'), 
Rt ={p^ + i?)i?^ 

+ e {PlRxxx + hRxxRx + PsRx) 

+ e (5li?5a; + S2RixRx + ^sRxxxRxx + S^R^^Rx + d^RxxxRx + ^qRxxRx + 1^7 -^a;) 

where, again, coefficients are certain functions of i? (it can be demonstrated that all odd order 
corrections in e must vanish identically). Substituting Eqs. (fT3]) into Eqs. ([TO]) , using Eqs. (dH) 
and the compatibility condition Ryt = Rty, we arrive at the recursive formulae for higher order 
corrections which are uniquely expressed in terms oi fi{R) and its derivatives. Thus, at the order 
e^ we obtain 

bi = ^, b2 = ^^{p"-{p'n 

^1 = ^, Pi = ^{2pfi' + l), 

cp, = 1(2^" - if,'f), P2 = ^iW - Mp'f + ^(p'f), 

^3 = Y^ip'" - Kp'fp"), /53 = ^{pp'" + 2/xV" - Mp'fp" - {p'f), 

which implies ([7]). Similarly, at the order e^ we have 

di = ^^{2,p" - {p'fl 

d2 = ^^{Qp"'-l%p'fp" + {p'f), 

d, = r^{^p"'-np'fp" + A{^,')% 

d, = 1^(11/^^'''^ - 49(/.') V" - W{p"f + 52iprp"), 

d, = ^(5^^"'^ - 32(/.')V('^) - 146mW" + 51(^0 V" - 44(/)3 + 216(^')'(/)'). 
The remaining coefficients can be expressed as follows: 

Pi = di, p2 = d[ + d2, P3 = d2 + 2^3, 

P4 = d'^ + 2d4, P5 = d2 + d4, p6 = d'^^ + Ms, Pr = d'^, 



5i = 2M + Y^(^')', 

62 = 2M + d2) + 3^(720(ii - 2(^0' + 5/), 

53 = 2^(d2 + 2d3) + Y^ilUOdi - 3{f,'f + 8fi"), 

J4 = 2^(4 + 2d4) + ^(2160di + V) + ^ (432^2 + 864^3 - 22(^') V" + V + (^0^) 

55 = 2^(4 + ^4) + ^(480di + fi") + ^ (720^2 - 16(^0 V" + V" + (^')') , 

^6 = 2^(4 + 4^5) + 2^"(d3 + 3^2) + J^^"'(rfi + 2^") 

zoo 

+ -^ (1440^4 + 18(m')V" - 78;u'(/i")' - 35(/i') V" + H^^'^^ 



57 = 2^4 + ^"d4 + ^ (di + 2/.") + ^(576(i2 + /i'") 

+ ^ f2304d5 + 21{fi'f{fi"f + 6(/x') V" - 48mV'V'" " 8(/x')V^'''^ " 24(^")' + 2^(^)) . 
576 V / 

Although we have calculated dispersive corrections up to the order e^ , the expressions are getting 
increasingly more complicated. It is important to emphasize that the calculation of dispersive 
corrections is an entirely algebraic procedure which does not require solving differential equations. 

2.2 Two-component reductions 

Two-component reductions of the dKP equation (jlip are sought in the form u = u{R^,R'^), 
w = w{R^,R'^), where the Riemann invariants R^ and R^ satisfy the two-component version of 
Eqs. ([3]). One has A* = ()u')^ + n, djW = ^WjU, where ^^ and u satisfy the two-component 
Gibbons-Tsarev system ^, 

r, 1 d2U Q 2 ^1'" o 2 diu d2U 

^^^=^73^' ^'^=]I^^^ ^''""^W^T?- ^^ 

We consider (1 + 1)— dimensional decompositions of the KP equation (|1U|) obtained by deforming 
Eqs. dSD, 

2( i ■ i ■ k i ■ k l\ 4 ^'^^^ 

R^ = X R^ + e iAjR^^^ + B^f^R^^R^ + C^j^iR^R^R^j + 0(e ), 

where the coefficients a^pA^-, etc, are certain functions of R^,R?] here i,j,k,l = 1,2. We point 
out that u{R^,R'^) and w{R^,R'^) remain undeformed (this can always be achieved modulo the 
Miura group). Moreover, one can show that all terms at the odd powers of e must vanish 
identically. Substituting u{R^,R^) and w{R^,R'^) into Eqs. p^ . using Eqs. (fT6|) and imposing 
the compatibility condition i?L = i?L, one uniquely recovers expressions for the coefficients in 
terms of /x* and u. Although these expressions are quite complicated in general, some of them 
are remarkably simple. Thus, 

1 1 difi^ 1 1 92Ai^ 2 1 ^lA*^ 2 1 ^2At^ 



^1 - T7:~a~T' ^2 - -TTT^r— ' oi - -7^-5--' ^2 



12 aiu ' ^ 12 5iu ' ^ 12 52U ' ^ 12 52U 



or, in tensor notation, a*- = j^i — lY'^^dj^^/diU. Similarly, 



A\ 



12 6diu ' ■ 
1 



ifl^ + is')d2fi' 



Al 



(/xl+ ^2)5^^1 



Ai 



12 652M ^ 



125iii ' 1292U 

or 4. = l^{{-lY+^{^ji + fi^)djfi^/d,u + 6j). 

As a particular case, let us consider deformations p6|) such that the series truncates at the 
order e'^. This leads to the following two possibilities: 

Case 1. 

diu= {fi"^ - fj.^)din^, d2U = {fi^ - fj.'^)d2fJ.'^ . (17) 

Eqs. (J17p together with the Gibbons-Tsarev system (J15p imply 

/i + /i = const. 

Setting fi^ = a + p, fi^ = a — p where a = const, and substituting this representation into the 
first two equations (jlSp . one obtains u = b — p"^, b = const. Then Eq. (jlSP '^ implies 9i2(p^) = 0, 
so that, up to a reparametrization of Riemann invariants, one has p = {R^ — B?Y'^ . Ultimately, 



+ {R'-R^f/\ 



a+{R^-R^fl\ 



fji = a + {w — w }''" , fj,' 

Up to a linear transformation of the independent variables, one can set a = b = 0. This implies 

A' = and R] = 0, so that the system (fT6]l 9 becomes trivial, while the system (fT6l) i reduces to 

2 
the stationary Boussinesq reduction of the dKP equation, {uux + f2''^a;xa;)a; + Uyy = 0. 

Case 2. 



diu = -(//^ -^2)(9i^\ 



1 



^2^^ = -(/U^ - /i^)92/x2. 



(18) 



Eqs. (118]) together with the Gibbons-Tsarev system (115]) imply that, up to a reparametrization 
of Riemann invariants, one can set 

1 



fJ' 



4 4 



/^ 



2-3^2^1^1 



4 4' 16 

The corresponding e^— coefficients in Eqs. (J16p take the form 

12 1 12 1 

«i=«i = TT777^ 75tY> 02 = ^2 



(i?i - R^) 



.2\2 



^1 



5i?i + i?2 



^, 



2(i?i-i?2)' 

1_ i?l+i?2 



Al 



2(i?2_i?i)' 



^2 



R^ + 5i?2 



6(i?i-i?2)' ^2 2{R^-R^y ^^' 2(i?i-i?2)' ^^2 6(^2 _^i)' 

etc. Introducing the new dependent variables v = {R^ + R^)/2 and u = [R^ — R'^)'^ /16, one can 
rewrite Eqs. (fT6]) as 

-Uy = {uv)x, 

2m2 



VVx + Ux + — 



Uxx 
U 



and 



ut 



Vt 



g2 / ^2 

(2n + f2)nx. + 2n7;'t;j. + — ( Au^x - 3^ 

12 \ u 



2vux + (2n + v'^)vx + 



n 



1 u2\ 1 / n^. 

o^ +^ 2v^^ + 3— V, 



respectively. This is the well-known Zakharov reduction of KP to NLS |32j . 



2.3 Waterbag reduction 

It was observed in |25j that n-component reductions of dKP can also be sought in the form 

u = u[v^,v'^, . . . , v"'), w = w{v^ ,v'^, . . . , f"), where the fields w* satisfy the equations 

4=(^ + ^) ' vl=(^^ + uv^ + w^ . (19) 

The system (I19p is automatically compatible provided u and w solve Eqs. (jlip . The substitution 
into Eq. (fTT]) 9 implies 

diW = V^ dkU + w* 5jn, (20) 

here di = d^i. It turns out that Eq. (fTT]) i is satisfied identically modulo Eq. ([201) . The 
consistency conditions didjW = djdiW imply [25] 



diU — djU + v^ — v^ 



9^3U = — ^ ^. , ^., ^., , (21) 



which can be viewed as an analogue of the Gibbons-Tsarev system ([5]). In the two-component 
case, n = 2, the system (PU]) simplifies to 

diw = {diu + d2U + v^)diu, d2W = {diu + d2U + v^)d2U, (22) 

and the compatibility condition did2W = d2diw takes the form 

d2U duu - diu d22U , . 

^12-" = 7^ t; — — ^ 2"; (23) 

OiU — O2U + v'- — v^ 

this is an analogue of Eqs. p^. The special case 

1 2 ai(t;-^)^ + a2(w^)^ , nz 1 2\ / .n 

n = aiu + a2V , w = h (ai + a2)(ait' +02^), (24) 

where ai and 02 are constants, is known as the waterbag reduction |17j . We seek a deformation 
of the two-component waterbag reduction p9|) in the form 

_ fiV? 






^^3 X (25) 



i = 1,2. Modulo the Miura group, we assume that Eqs. ()24p remain undeformed. Substituting 
Eqs. (p^ and (p5]) into ([TT]) and using the compatibility conditions Vy^ = vly, we obtain explicit 
expressions for P* and Q^. Thus, Eqs. ()25|) i take the form 



yl=l^J^ + u] +£^P + 0(e4), 



/fr;2)2 X ,2 (26) 



" \ 2 L, a2 



where 

^ A + 2ai - 02 1 A + ai - 2a2 o 

^ 1 9 A XXX 1 9 A a;a;a; 

m — a2A(A — ai — 2a2) i i nA + m i n 2 i\ m + aiA(A + 2ai + 02) 2 
^ 12aia2A2 ^^"^^ 12aia2A2 ^^^'^^ ^ '^^^'^^^ ^ 12aia2A2 ^^^ 

24aia2A3^ ^^ 24aia2A3 ^ ^^ ^ 24aia2A3 ^^ ""^ 24aia2A3^ ""^ 
here 

A = V — V , 171 = ai{ai — 202) — a2('^2 — 2ai), n = a^^ + 0102 + a2- 

Notice that, although the hydrodynamic part of Eqs. (j26|) is conservative, the expression P is 
not a total x-derivative. The expressions for Q^ are not presented here due to their complexity. 

3 Deformations of the dTl reductions 

Expanding the r.h.s. in Eqs. ([8|) one obtains 



Uy/u = Wx- -Wxx + -prWxxx H h (-l)""*"^— rtfnx + 

" 2 6 n\ 



e e2 £■ 



(27) 



li^i — Ux + ri ^xx ~r Uxxx + ' ' ' + . "U-na:; + 

2 6 n! 



The corresponding dispersionless limit, the dTl equation, is 

Uy = uwx, wt = Ux. (28) 

It admits one-component reductions of the form 

u = R, w = w{R), 
where w' = l//x, and R{x, y, t) satisfies a pair of Hopf-type equations 

Rt = fJ-Rx, Ry = — Rx] (29) 

here n{R) is an arbitrary function. Deformations are sought in the form 

u =R, 

w =w{R) + eaRx + e^ihRxx + 62 -R^) + e^iciRxxx + CiRxxRx + C3RI) + O(e^), 

where a, bi, Ci are certain functions of R, and 

Rt =^lRx + e\piRxxx + P2RxxRx + PsRl) + O(e'), 

R (31) 

Ry =-Rx + e^i^iRxxx + f2RxxRx + fsRl) + O(e^). 



(30) 



We point out that the relation u = R remains undeformed (modulo the Miura group), and that 

all odd order corrections in e in the equations (fSTI) vanish identically. Substituting (pO|) into (f27|l , 

using (j3ip and the compatibility condition Ryt = Rty, one can recursively calculate all higher 

order corrections in terms of fi and its derivatives at different powers of e: 

Order e: 

1 

a = — ; 
2//' 

Order e^ : 

Ru' - 2/i 
oi =- 



b 



'2 



12^2 ' 



24;u4 

¥^2 =Y^ (2i?^(^')' + ^'(^')' - m'(3^' + 2i?/i")) , 

</^3 = - ^ {4R\f^T - 2mV'(«m" + mO - Rfi{f^')\3Rfi" + 2/x') + fi^Rf,'" + 2/)) ; 

Order e^: 

Rfj,' — fi 



^i = -titt:2(^^'-^)' 



Cl 



24/i2 ' 



C3 = - ^ (4ii2(/u')' - 2^?^^'{R^^" + ^') - Rf^ifi'fiSRfi" + 2/i') + /.^(i?^"' + 2;,")) , 
etc. Again, this calculation is an entirely algebraic procedure. 

4 Recostruction of dispersive terms in 2D 

Given a dispersionless integrable system in 2 + 1 dimensions, how can one reconstruct the 
corresponding dispersive counterpart? This natural problem was first addressed by Zakharov in 
[33], based on the quantization of the corresponding dispersionless Lax pair. This, however, does 
not work when the Lax pair is 'sufficiently complicated'. We propose an alternative approach to 
this problem, based on the requirement that all hydrodynamic reductions of the dispersionless 
system should be inherited by its dispersive counterpart. This imposes strong restrictions on 
the structure of dispersive terms. The following simple examples illustrate the method; further 
results in this direction will be reported elsewhere. 

10 



4.1 Generalized KP equation 

A two-component form of the generalized KP equation ([9|) is 

ut - uux + e{AiUxx + Mul,) + e^{BiUxxx + B2UxUxx + Bsul) = Wy, 



(32) 



where Ai and Bi are certain functions of u. Let us require that any one-component reduction 
([6]) of its dispersionless hmit can be deformed, in the form (J14p . into a reduction of Eqs. (j32p . 
Thus, whenever one obtains a relation involving fi{R) and its derivatives, all coefficients must 
be set equal to zero: this ensures that n{R) remains arbitrary. Looking at different powers of e 
we arrive at the following results: 
order e : 

Ai = 0, A2 = 0; 

order e^: 

B' 
Bi = const, S3 = — ; 

order e^: 

WiB'^ - 2Bl = 0. 

This implies either B2 = 0, or B2 = —3Bi/{2u + cq), where cq is an arbitrary integration 
constant. It follows that only the case 

^2 = 0. 



is admissible. Thus, the KP equation is the only possible extension, within the class ()32p . which 
inherits deformations of all hydrodynamic reductions. 

4.2 BKP and CKP equations 

Let us consider equations of the form 

Ut - 5(n^ + w)ux - 5uwx + 5'Wy + e^{AQUUxxx + AiWxxx + A2Uxxx) + e^^s Uxxxxx = 0, 

{66) 

Uy = Wx, 

where Aq — A^ are arbitrary constants. Deformed one-component reductions are sought in the 
form 

u =R, 

w =w(R) 

9. 9. (34) 

+ e\hRxx + b2Rl) 

+ e {diR^x + d2RxxxRx + d?,Rxx + d^RxxRx + d^Rx) + 0{e ), 
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(35) 



where R obeys the system of the form 

y — l-^^x 

+ e {'■PlRxxx + ^2RxxRx + ^sRx) 

+ e {PlRbx + P2RaxRx + PsRxxxRxx + PaRxx^x + P^RxxxRx + PqRxxRx + P7-Rx] 

+ 0(e'), 
i?t=5(i?2 + u> + ^(i?-;u))i?^. 

+ e (/?1-Ra;a;a; + P2RxxRx + Pz^x) 

+ e (^l-Rso; + 52RaxRx + ^sRxxxRxx + ^aRxx^x + ^bRxxxRx + ^qRxxRx + ^7-^0;) 

+ 0(e'); 

here /^(i?) is an arbitrary function, and if' = /i. Substituting Eqs. (fM|) and (f35]) into Eqs. ([33]) . 
we obtain (/?, p, f3, 5 in terms of b, d, etc. The compatibihty condition Rty = Ryt provides the 
following constraints on the constants j4j's, 

order e^ : 

A^ = An, 



order e^ 



A^ = --S, ^Al - TAnAo + 2Al = Q. 



25 



The quadratic equation leads to either A2 = Aq, or A2 = \Aq, which correspond the BKP and 
CKP equations, respectively. We have verified that, up to the order e^, no extra constraints 
appear, and that all coefficients in the expansions ()34p and (I35p are certain expressions in terms 
of p. This calculation shows that the reconstruction of dispersive corrections is essentially non- 
unique. 

4.3 The 'universal' equation 

In this section we study dispersive deformations of the so-called 'universal' equation, 

(36) 

Ut =Wy + UWx — WUx, 

which was discussed in a series of publications [24^ l26l |3], etc. The structure of hydrodynamic 
reductions in this case is particularly simple, namely, n-phase solutions are given by the formulae 

u = R^ + ... + i?", w = f\R^) + ... + r(i?") - ^{R^ + ... + i?")2, 

where f^{R^) are n arbitrary functions of one variable, and the phases -R* satisfy Eqs. ([!]) with 

notice that ^* = diW. In the one-component case one has u = R, w = w{R), where the single 
phase R satisfies the equations 

Ry = p{R)Rx, Rt = {p'^ + Rp-w)Rx, (37) 

12 



w' = //. We have found the following dispersive deformation of Eq. 

Uy =Wx, 

4 8 (38) 

Ut =Wy + UWx - WUx + re {UxxWxxx - WxxUxxx) + 0{e ), 

which inherits deformations of all one-components reductions (|37p up to the order e^. Our 
procedure, however, does not work beyond this order, leading to the trivial condition r = 0. We 
would like to emphasize that the extension of two-component reductions leads to the condition 
r = already at the order e^. The deformation (f38|) is found by seeking dispersive extensions 
of Eq. (j36p in the form 



^y — ^xi 

Ut =Wy + UWx — WUx 

+ Ae{uwxx - wuxx) 

+ Bie^{uWxxx - WUxxx) + B-2e^{UxWxx - WxUxx) (39) 

+ Cie^{uWix - WUix) + C2e^{UxWxxx - WxUxxx) 

+ Die^iuW^x - UJU5x) + D2e^{UxWix - WxUAx) + D'ie^{UxxWxxx - WxxUxxx) 

+ Eie^{uWQx - WUQx) + E2e^{UxW5x - WxU^x) + Ese^iUxxWix - WxxUAx) + • • • , 

where the coefficients A, Bi, etc. are constants. Notice that one can simplify Eq. (|39p by using 
Miura-type transformations of the form 

u —> u + aeux + /?e Uxx + •••) w ^ w + aewx + /3e Wxx + •••) 

where a, (3, ... are arbitrary constants. In particular, using this freedom one can set B2 = C2 = 
D2 = E2 = 0, etc. In what follows, we assume this normalization. Let us require that any 
one-component reduction (1371) can be extended as 

u =R, 

w =w{R) 

+ eaRx + e^{biRxx + b2Rl) (40) 

+ e {ciRxxx + C2RxRxx + C'iRx) 

+ e (diRix + d2RzxRx + dsRxx + d^RxxRx + d^Rx) + • ■ ■ , 
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(41) 



where R satisfies the system 

Ry =nRx 

+ e{aiRxx + <T2-R^) 

+ e^{ipiRxxx + ^2RxxRx + VsRI) 

+ e^^lpiRix + tp2RxRxxx + ipsRlx + ipiRxxRl + tp^Rt) 

+ e [piRbx + P2RixRx + PsR-SxRxx + PaRxxR-x + P^R^xRx + P&RxxRx + PlRx) 

+ ..., 

Rt ={p. + Rp — w)Rx 

+ e(aii?a;a; + a2-R^) 

+ e {PlRxxx + hRxxRx + PsRx) 

+ e^(7ii?4x + l2RxRxxx + 73-RL + liRxxRl + 75^x) + ^aRxxRI + V'S^x) 

+ e ((5ii?5a; + 52RixRx + (^S-Rsx-Rxa; + ^^RxxRx + ^bRsxRx + ^qRxxRx + <^7-Ra:) 

+ ..., 

tf' = /x. Substituting (PUJ) into ([5U]) . using Eqs. ([^T]) and the compatibihty conditions i?yi = i?ty 
(where p is assumed to be arbitrary), at different orders in e we get the following conditions for 
the coefficients of ([5^ . 
order e: 

order e^: 

order e^: 

order e^: 

order e^: 



A = 0, 

Ci = 0, 

El = E'i = 0, 



etc. Notice that there is no condition on the coefficient D^ = t. This results in Eq. (|38p . The 
corresponding Eqs. (I40|) . (I4ip take the form 

w =w{R) 

+ e [diRix + d2R3xRx + dsR^j-^j. + d4RxxRx + d^Rx) + 0(e ), 

and 
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Ry —flRx 

+ e {piR^x + PlRixRx + PsRsxRxx + PaRxxRx + PbRsxRx + PqRxxRx + PtRx) 

Rt ={p^ + Rp- w)Rx 

+ e ((^l-Rsx + ^2RixRx + ^sRsxRxx + 5iRxxRx + S^RsxRx + ^qRxxRx + <^7-Rx) 

respectively. Here all coefficients are certain explicit expressions in terms of /i and its derivatives, 

r^-- 3r/.-- t(3/." + 4)/.'" 

di = 0, d2 = TT-;; — ^, as = -tt";; — ^, "4 " 



3;u"-2' ^ 3;u"-2' 3(;u")2-llAi" + 6' 

_ r ((74 - IV - 12(/i")2)(^'")^ + (6 - 23/u" - 17(/i")^ + S{p"f)p'^^^y) 
'^^~ 8(/i"-3)2(2-5^" + 3(^")^) ' 

(Jl = 0, (^2 = (i? + 2/i')p2, ^3 = (i? + 2;U')P3, 

5i = [R + 2^')P4 + 3^"p3 + 3r^" - ds, 

^5 = (ii + 2^')P5 + ^P"d2 - Tp" - d2, 

Se = {R + 2/i')p6 + 2/i"'(3d2 + ds) + r/^'" + 5p"d4 - ^4, 
Jr = (i? + 2/x')/07 + P^^^'^d2 + /x"'(i4 + ^p"d5 - 4, 

/Ol =0, P2 = ^2, P3 = C?2 + 3(i3, 

P4 = (ifj + 2d4, P5 = d2 + d4, pq = d'^ + Ad^, Pt = d'^. 

Remark. The paper [3] provides a multi-parameter deformation of the universal hierarchy with 
the first few terms given by 



€ LL € U 

+ ^Tg- (Wa'xxW^4x - 'Wxxx'U4a;) " "^Tg" WxxW^x - WxxU^x) + • • • , 



Notice that this deformation belongs to the class (j39|l . According to our calculations, it does 
not inherit hydrodynamic reductions of the universal equation. 

5 Concluding remarks 

Several natural questions arise in the present context: 

1. The results of (HO [21] suggest that any evolution system in 1 + 1 dimensions is quasitrivial, 
that is, transformable to its dispersionless limit by a quasi-Miura transformation. In particular, 
deformed hydrodynamic reductions constructed in this paper are quasitrivial as well. It would 
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be of interest to extend these quasitriviality type results to the corresponding integrable systems 
in 2 + 1 dimensions such as the full KP and Toda lattice equations. 

2. As demonstrated in [33], dispersionless limits of integrable systems in 2 + 1 dimensions come 
with Lax pairs of the form 

7py = A{u,^P^), Vt = ^(u,Vx), (42) 

which generate equations for the fields u(x, y, t) via the compatibilty conditions 

At - By + {A, B} = 0; 

here {^4, B} = A^B^ — B^A^ denotes the standard Poisson bracket of A and B, and ^ = ipx- It 
was shown recently (see e.g. [9tllH[T3]) that, conversely, for broad classes of multi-dimensional 
dispersionless systems, the existence of Lax pairs of the form (I42p is necessary and sufficient for 
the integrability. In general, the dependence of A and B on ^ can be rather non-trivial, which 
makes a direct 'quantization' of such Lax pairs difficult. Several ways to solve the problem 
of quantization of dispersionless Lax pairs of the form ()42p were proposed in the literature, in 
particular, including the Moyal deformation of the Poisson bracket { , }, see e.g. [331 [2S1 [2H1 El H] • 
These methods work well for a limited number of examples, allowing one to reconstruct the 
admissible dispersive terms, however, they meet difficulties when applied to Lax pairs with a 
more complicated dependence on S^. Our approach to the problem of quantization is to allow only 
those dispersive corrections for which all hydrodynamic reductions of the dispersionless system 
can be deformed into reductions of its dispersive counterpart. We have demonstrated that this 
requirement is very restrictive indeed, and imposes strong constraints on the structure of possible 
dispersive terms. It would be challenging to apply this recipe to other classes of dispersionless 
integrable systems obtained in the literature, and to classify the associated soliton systems in 
2 + 1 dimensions. 

3. Our approach provides an infinity of decompositions of a given (2 + l)-dimensional integrable 
soliton equation into a pair of commuting (1 + l)-dimensional flows, which are parametrized 
by arbitrarily many functions of one variable. An alternative construction is known as the 
method of symmetry constraints, or potential-eigenfunction constraints, see e.g. |18 1I191 [20] and 
references therein. It would be of interest to understand whether generic deformed hydrodynamic 
reductions possess any kind of 'symmetry' interpretation. 

4. Solutions of dispersionless systems coming from hydrodynamic reductions are known to break 
down in finite time. The addition of dispersive corrections can be viewed as a regularization 
procedure preventing the gradient catastrophe. Although the structure of higher order correc- 
tions can be rather complicated in general, it is the very first term that seems to be of prime 
importance. Thus, our procedure gives a canonical way to regularize hydrodynamic reductions. 
It would be of interest to investigate the behavior of regularized solutions numerically. 
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